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In this work we examine the Cauchy convergence of both post-Newtonian (T-approximant) and re- 
summed post-Newtonian (P-approximant) templates for the case of a test-mass orbiting a Kerr black 
hole along a circular equatorial orbit. The Cauchy criterion demands that the inner product between 
the n and n + 1 order approximation approaches unity, as we increase the order of approximation. 
In previous works, it has been shown that we achieve greater fitting factors and better parameter 
estimation using the P-approximant templates for both Schwarzschild and Kerr black holes. In this 
work, we show that the P-approximant templates also display a faster Cauchy convergence making 
. them a superior template to the standard post-Newtonian templates. 
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I. INTRODUCTION 

Or 

The emission of gravitational waves from double neutron stars (NS), double black holes (BH) or mixed binaries of a 
neutron star and a black hole, are expected to be major detectable sources for the interferometric detectors, LIGO 
VIRGO @, GEO600 0, TAMA and LISA Due to radiation reaction, the binary orbit slowly decays emitting 
| a "chirp" signal (i.e. a signal whose am plit ude and frequency increases with time). While there is probably a greater 
population of NS-NS binaries 0, Q, H, it is BH-BH binaries that are the strongest candidates for detection since 
they can be seen from a greater volume |o. Ill|. In order to detect such sources one employs the method of matched 
filtering |12j . In this method, a set of theoretical waveforms or templates, that depend on a number of parameters 
of the source and its location and orientation relative to the detector, are cross-correlated with the detector output 
weighted by the inverse of the noise spectral density. The success of matched filtering depends on how well the phase 
■ evolution of the waveform is known. To improve the recovered SNR for inspiral events, there have been many efforts 
to accurately compute the dynamics of a compact binary and the waveform it emits, or, to use phenomenologically 
defined detection template families [l3L IlH ITo| . 

Two different approximation schemes have been used to try and accurately predict the phase of the GW. On 
one hand, there is the post-Newtonian (PN) expansion of Einstein's equations [lti, liH ITsl IT9I I20L l2ll Izl which 
assumes slow motion bodies. On the other hand, there is black hole perturbation theory |23l l24l l25l l26l \2li. l28j 
. . , which has no restriction on velocity, but assumes that one body is much more massive than the other. Due to the 
fact that the convergence of both post-Newtonian ap pro ximation and black hole perturbation theory seem to be too 
slow to be useful in constructing accurate templates |28l l29l 1.3(1. 1.'i ll. |32| . Damour, Iyer and Sathyaprakash (hereafter 
DIS) showed for the case of a test-mass in orbit about a Schwarzschild BH, that newly constructed P-approximant 
templates achieved larger overlaps with the exact waveform and lower biases in the estimation of parameters than 
the corresponding post-Newtonian (hereafter T-approximant) templates. While more P-approximant templates are 
needed to cover the same volume of parameter space (33], the increased performance of the template warrants the 
extra computational cost. Most recently, the DIS analysis was extended to the case of a test mass orbiting a Kerr 
BH (34|. It was shown, once again, that the P-approximants outperformed the T-approximant templates. For a more 
indepth discussion on the definition of the gravitational waveforms and the construction of the P-approximants, we 
refer the reader to Reference l34l. 
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II. THE CAUCHY CONVERGENCE OF T- AND P-APPROXIMANT TEMPLATES. 

If a sequence {x n } converges, the terms get closer and closer to the limit of the sequence as the order of the 
approximation increases. The Cauchy Convergence Criterion states that if and only if for every e > 0, there exists 
an integer N , such that for each pair of integers m > N and n > N , \x m — x n \ < e. The main difference between the 
Cauchy criterion and other convergence criteria, is that instead of demanding that the terms get closer to some limit, 
we demand that the terms get closer to each other. We therefore require that for our approximations to be Cauchy 
convergent, as we increase the order of PN approximation we demand the condition 

(h n \h n+ i) — ► 1 as n — > 00, (1) 

where 

^) = 2 J o °°s^f) [HfW(f) + h*V)9(f)\, (2) 
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the * denotes complex conjugate and h(f), g(f) are the normalized Fourier transforms of h(t), g(t). We can see that 
the inverse of the detector one-sided power spectral density, Sh{f), acts as a weight in the definition of the scalar 
product. In the case of circular equatorial orbits of a test mass around a central black hole, the wave is parameterized 
by the time of arrival of the wave to, the phase at time of arrival 0o, the total mass m, the reduced mass ratio 77 and 
the spin of the central black hole q. In calculating the Cauchy convergence, the masses and spins of both templates 
are kept the same. We therefore only need to maximize over the parameters to and <po. Maximization over to is 
achieved by simply computing the correlation of the template with the data in the frequency domain followed by the 
inverse Fourier transform. This yields the correlation of the signal with the data for all time-lags. 

Maximization over (j)Q is achieved by what DIS define as the " Best Possible Overlap" when individually maximizing 
over the phases of two separate templates : 



(3) 



where 
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and 



where = h((j)o = 0) and hf 
possible overlap in all cases. 

In this paper, as we are working in the test-mass approximation, we assume that our system is composed of objects 
with a small mass ratio. We have thus chosen the following five systems to analyse : (100,1.4)M Q , (50,1.4)M Q , 
(10, 1.4)M , (100, 10)M Q and (50, 10)M Q . While our results are more reliable for the more extreme mass ratio cases, 
it is also in our interest to examine the behaviour of templates as we approach a quasi-equal-mass case. It was 
shown in reference |34| that singularities plague the sub-diagonal P-approximant at the x 7 order, i.e. P% , making this 
template completely useless. Such singularities do not exist at the same order for super-diagonal P-approximants, 
i.e. P3 . Because of this, when we examine the Cauchy convergence involving the x 1 approximation we will use the 
super-diagonal P-approximant as our template. In this work we will use the projected third-generation EURO noise 
curve in our overlap calculations as it has the greatest sensitivity |35j . 

III. RESULTS AND DISCUSSION. 
A. Black Hole - Neutron Star. 

In Figure (JIJ we have plotted the Cauchy convergence for T and P-approximants for the systems (100-1.4), (50-1.4) 
and (10-1.4) M Q . In the high mass case we can see that the T-approximants have excellent convergence for the 
extreme retrograde spin at q = —0.95. Once we increase the spin value and move towards prograde orbits we can 
see that these templates do not display good convergence properties, with the extreme prograde case at q = 0.95 
never giving an overlap of > 0.2. As we move down the column on the left hand side of Figure Q we can see that 
the Cauchy convergence gets progressively worse as we move towards the equal mass case. In all cases we find that 
the best convergence is achieved for retrograde orbits. This should not surprise us as these are slow motion systems. 
Thus, the overlaps get expectedly worse as we increase spin and decrease the total mass of the system. In fact, in the 
case of a (10-1. 4)M system, we can see that besides a spin of q = —0.5, we obtain very bad convergence at all spin 
levels with the T-approximant templates. 

In contrast, the P-approximant templates in the (100-1. 4)M Q case show excellent convergence properties. Again 
we can see that the convergence begins to tail off as we increase in spin. As we move down the cells on the right 
hand side of Figure (Q) toward the equal mass case, we can see that the P-approximants have slow convergence at 
low orders of approximation, but approach unity as we move towards the a; 8 approximation. This again is expected 
behaviour as we know that the P-approximants do not perform as well at lower orders of approximation. We can 
however see that even in the (10-1. 4)M Q case we achieve overlaps of > 0.9 up to a spin of q = 0.5. We should remark 



(cos 9ab) t 



A + B 




A-B 



C 2 



= (et\e?)+(e?\e*), B = (ef |ef > + (ef |ef > , (4) 
= <ef|ef) <4|ef>+(ef| e f) (e£ |ef ) , 

a R_ p a = H-W\<*) <± ~ 
1 1^1 ' 2 \h£-(h£\et)e?\ (5) 

= h((f>o = ir/2). When we calculate the Cauchy convergence, we will quote the best 
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FIG. 1: The Cauchy convergence of the T and P-approximant BH-NS systems (100-1.4), (50-1.4) and (10-1.4) M for the 
spin values q — —0.95 (circle), q = —0.5 (square), q = (diamond), q — 0.5 (triangle), q = 0.95 (upturned triangle). The 
T-approximants occupy the column on the left, while the corresponding P-approximants are on the right. All value shown are 
the Best-Possible-Overlaps between templates. 




FIG. 2: The Cauchy convergence of the T and P-approximant BH-BH systems (100-10) and (50-10)^/© for the spin values 
q — —0.95 (circle), q = —0.5 (square), q = (diamond), q = 0.5 (triangle), q = 0.95 (upturned triangle). The T-approximants 
occupy the column on the left, while the corresponding P-approximants are on the right. All value shown are the Best-Possible- 
Overlaps between templates. 



here about the convergence properties of both templates for the extreme prograde case. Graphically, it looks like 
the P-approximants perform only marginally better than the T-approximants. However, we mentioned before that 
the flux function goes to zero for the T-approximants before the LSO is reached. Thus, these templates are shorter 
than they should be and never reach the truly relativistic regime of the waveform. The P-approximants on the other 
hand have been evolved to the LSO, so the convergence results are more reliable. We cannot have confidence in the 
results from the T-approximants in the extreme prograde case, because if we could evolve these templates to the LSO 
we would expect the convergence to be even slower. The P-approximants thus produce higher overlaps with a full 
template which point towards their superiority as templates. 
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B. Black Hole - Black Hole. 

In Figure we have plotted the convergence results for the systems (100-10) and (50-10)M©. If we again focus 
on the cells in the left hand column, we can see that we do not have the same pattern of performance for the T- 
approximant templates. For the (lOO-lO)M© case, we can see the T-approximants obtain the fastest convergence as 
we move towards the Schwarzschild case. Once we increase the spin of the central BH and move onto prograde orbits, 
we again see a drop-off in the convergence rate. In the (50-10)Af Q case we can see that the Cauchy convergence 
gets progressively worse as we move from retrograde to prograde motion. The best overlaps are again obtained from 
extreme retrograde orbits, while the worst are from the extreme prograde. 

Once again, focusing on the cells in the right hand column, we can see that the P-approximants have better 
convergence properties. While in the retrograde case for the (lOO-lO)M© the convergence is slower, we can see 
that for both systems the overlaps approach unity as we move to higher approximation. In most cases, while the 
convergence is slower at lower orders of approximation, it is still accelerated in comparison to the T-approximant 
templates. Again we can see that while displaying a faster convergence, the P-approximants at q = 0.95 still have a 
poor Cauchy convergence. While it is obvious that the P-approximants are superior templates to search for waveforms 
up to a spin of q ~ 0.5 (and possibly slightly higher), they may still not be the best templates to search for extreme 
prograde systems. 

IV. CONCLUSIONS 

We have looked at the Cauchy convergence of both T and P-approximant templates in the test-mass approximation 
for Kerr binaries. We examined the BH-NS and BH-BH cases for both prograde and retrograde orbits using the third 
generation detector, EURO. In a previous work it was shown that P-approximant templates achieve higher fitting 
factors with lower biases in the estimation of parameters compared to the T-approximants. In this study we have 
shown that in all examined cases with various spin values, the P-approximants have a faster Cauchy convergence than 
the T-approximant templates. We found that the T-approximants show good convergence for extreme mass-ratio 
retrograde motion. However in all cases there is a serious retardation in convergence as we move to Schwarzschild and 
prograde orbits. We also see a degradation in overlaps as we move towards the equal mass case. The P-approximants, 
while slow to converge at low orders of approximation, approach unity in all cases as we increase the approximation 
order. This is due to the known fact that the P-approximants perform better at higher orders of approximation. 

While convergence is slow in both cases for the extreme prograde case, we have more faith in the P-approximants 
due to the fact that we can evolve these templates to the LSO. This isn't possible with the T-approximants as the 
flux function goes to zero before the LSO is reached. In fact, these templates should not be used to search for systems 
with a spin of q > 0.5. This work reinforces the conclusion of an earlier work that P-approximant templates are more 
reliable than T-approximants. They give better fitting factors, have lower biases in the estimation of parameters, and 
as shown in this work, display an accelerated Cauchy convergence. Due to these results, these templates should be 
used in any data-analysis strategy to detect gravitational waves from binary systems. 
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